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A magnetic field applied to a three-dimensional antiferromagnetic metal can destroy the long- 
range order and thereby induce a quantum critical point. Such field-induced quantum critical 
behavior is the focus of many recent experiments. We investigate theoretically the quantum critical 
behavior of clean antiferromagnetic metals subject to a static, spatially uniform external magnetic 
field. The external field does not only suppress (or induce in some systems) antiferromagnetism but 
also infiuences the dynamics of the order parameter by inducing spin precession. This leads to an 
exactly marginal correction to spin-fluctuation theory. We investigate how the interplay of precession 
and damping determines the specific heat, magnetization, magnetocaloric effect, susceptibility and 
scattering rates. We point out that the precession can change the sign of the leading x/T correction 
to the specific heat coefficient c{T) /T and can induce a characteristic maximum in c(T) /T for certain 
parameters. We argue that the susceptibility x = dM/dB is the thermodynamic quantity which 
shows the most significant change upon approaching the quantum critical point and which gives 
experimental access to the (dangerously irrelevant) spin-spin interactions. 



The study of quantum phase transitions is currently a 
very active field of research in theoretical as well as ex- 
perimental condensed matter physics. Particularly in a 
large number of metals - mostly heavy Fermion or tran- 
sition metal compounds - the critical fluctuations associ- 
ated with a quantum phase transition induce anomalous 
behavior in thermodynamic and transport quantities like 
diverging specific heat coefflcients or a linear resistivity 
quite distinct from the behavior of a conventional Fermi 
liquid. 

Experimentally there are three main methods to tune 
a system towards a quantum critical point: doping, pres- 
sure, and magnetic field. Doping has the disadvantage 
that it induces disorder and it cannot be easily adjusted 
within a single sample. These problems are absent if 
pressure is used as the control parameter of the quan- 
tum phase transitions. However, the presence of a pres- 
sure cell makes many experiments difficult. For this rea- 
son, many recent experiments^^^"'^ investigate field-tuned 
quantum critical behavior, where an external magnetic 
field is used to control the distance from the quantum 
critical point. Generally it is expected that the presence 
of a magnetic field changes the universality class of the 
transition as in its presence time reversal invariance is 
broken. In this paper, we will therefore analyze theoret- 
ically the quantum critical behavior of a clean itinerant 
antiferromagnet in three dimensions subject to a static, 
spatially uniform external magnetic field B. 

Such a situation has been investigated in a num- 
ber of experiments^"^. For example in CeCu5.2Ago.8^ 
and CeCu5.8Auo.2^ magnetic order can be suppressed by 
moderate magnetic fields. In these systems the quantum 
critical behavior induced by a magnetic field B appears to 
be qualitatively different compared to the critical proper- 
ties for vanishing field (controlled by pressure or doping) . 
In the presence of a field these systems seem to follow^ the 
predictions from spin-fiuctuation theory^^"^^ for three- 
dimensional nearly antiferromagnetic metals, while this 
is not the case for B ~ 0^^. Similarly experiments'^ in 



field tuned YbCus-^^Ala; appear to be consistent with 
spin-fluctuation theory, which is not found to be the case 
in YbRh2Si2 where magnetic order is suppressed by tiny 
magnetic fields'*. Recently, in CeCoIns^''' the supercon- 
ducting order was suppressed by a magnetic field - it is 
at the moment a controversial question whether the ob- 
served anomalous behavior is related to a superconduct- 
ing quantum critical point or whether magnetism plays 
a role in this system. 

Another interesting class of systems are insulators like 
TlCuCl3^ SrCu2(B03)2^ or BaCuSi206i° where antifer- 
romagnetic order has been induced by the application of 
a magnetic field B. These transitions^ can be interpreted 
as a Bosc-Einstein condensation (see below) of spin-1 ex- 
citations. The energy of the "spin-up" component of such 
triplets is lowered by B until it condenses at a critical 
field, B = Be, thereby inducing antiferromagnetic order 
perpendicular to the magnetic field. 

In contrast to classical transitions, the dynamics, 
i.e. the temporal quantum fiuctuations, of the order pa- 
rameter determines the nature and universality class of 
a quantum phase transition. For example, at the critical 
point of an insulating antiferromagnet, the dynamics of 
the order parameter $ can be described*^ as in a Klein 
Gordon equation {df — V^)*!". In such a system, typical 
frequencies lu scale linearly with the momentum, uj oc , 
where z = 1 is the dynamical critical exponent. In con- 
trast, in a metal the excitation of particle- hole pairs leads 
to a Landau damping*^ of the antiferromagnetic order 
parameter, {dt + V^)*!" and therefore z = 2. Here we 
assumed that the ordering vector Q is sufficiently small, 
Q < 2kp, such that low-energy particle-hole pairs with 
momentum Q exist. 

A magnetic field will have two main effects; first it 
will suppress (or in some cases^ also induce) magnetic 
order. More interesting is the second effect: It induces a 
precession of the magnetic moments S perpendicular to 
the magnetic field 

dtS^BxS (1) 
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FIG. 1: Schematic phase diagram of a quantum phase tran- 
sition with a control parameter r oc B — Be. 

and therefore modifies the dynamics of the order param- 
eter. The linear time-derivative also translates to a dy- 
namical exponent z = 2, and therefore the question arises 
how the precession competes with damping in a metal, 
which is characterized by the same z. For insulating sys- 
tems the physics of the precession term has been widely 
discussed The corresponding quantum critical be- 
havior as a function of magnetic field of such an insulating 
magnet in an external field is actually well known: it is 
expected to be in the same universality class as the quan- 
tum phase transition of a low-density interacting Bose 
Einstein condensate as a function of chemical potential. 
The linear time derivative idt^ of the Schrodinger equa- 
tion can in this case be identified with the precession 
term (1), see below. 

In the following, we study the interplay of ohmic damp- 
ing and spin precession terms in the case of a nearly an- 
tiferromagnetic metal. First we present the model for 
the order parameter field and a short derivation of the 
effective action. Then we list the renormalization group 
equations for the parameters of the model and use them 
to derive the behavior of the correlation length. In the 
following sections we calculate the specific heat, thermal 
expansion, magneto caloric effect, and susceptibility. We 
show for example that sufficiently large magnetic fields 
can induce sign changes in the critical contribution to 
the specific heat and that the susceptibility is particu- 
larly suited to probe the vicinity of the quantum critical 
point. Finally, we investigate the influence of the B-field 
on the scattering rate of the electrons. 

I. MODEL AND EFFECTIVE ACTION 

Following Hertz we describe the critical behavior of 

an antiferromagnetic metal entirely in terms of an effec- 
tive Ginzburg-Landau- Wilson theory of an order param- 
eter field ^{r,t) which represents the fluctuating (stag- 
gered) magnetization of the system. 

In the absence of a magnetic field, the quadratic part 



of the action takes the form^^ (assuming negligible spin- 
orbit coupling) 

^^[*] = ^/(0^**^'" + '' + l""l)*' 

where r measures the distance from the quantum critical 
point and momenta k are given relative to the antiferro- 
magnetic ordering wave vector Q. The |a;„| term arises 
from the (Landau-) damping of the spin-fluctuations by 

gapless fermionic excitations in the vicinity of points on 
the Fermi surface that are connected by Q (assuming 

Q < 2kF). 

How will this effective action change in the presence 
of a magnetic field? First, r = r{B) will acquire a mag- 
netic field dependence, for example, r will grow for larger 
fields in systems where antiferromagnetism is suppressed 
by B. Second, the magnetic field breaks the rotational 
invariancc and components of $ parallel and perpendic- 
ular to B will have different masses, and r_L, respec- 
tively. Third, as argued above, the magnetization will 
preccss around B; this is described by an extra term (in 
coordinate and time space for convenience) 



5f [*] = JdT Jdrb-ii^ X dt^) (3) 



in the effective action, where b is parallel to B (taken to 
point into the ^-direction) and we have introduced the 

complex field $^ = ^x + i^y. Note that (3) breaks time- 
reversal invariance. Therefore such a term is absent for 
B = 0. 

Above, we deduced the form of the effective action on 
phenomenological grounds but it can also be derived from 
a more explicit calculation starting from a Hubbard-type 
model of electrons moving in the presence of a magnetic 

field, H = Y.kai'^k + Ba^„)'tl)l^'tpka + UY^nin-^. Here, 
the magnetic field enters only via a Zeeman-term; we 
do not take orbital effects into account, assuming that 
Landau levels are broadened by disorder or thermal ef- 
fects. Note that in the experimentally most relevant 
heavy Fermion system, orbital effects are strongly sup- 
pressed compared to contributions from the Zeeman term 
as the effective masses and magnetic susceptibilities are 
very large in those systems^^. 

For simplicity, we assume that the antiferromagnet is 
commensurate (incommensurate antiferromagnets show 
the same qualitative behavior for all quantities discussed 
below) and introduce a real order parameter vector 
^{x,t) as a Hubbard-Stratonovich field which decou- 
ples the spin-density part of the interaction. Follow- 
ing Hertz, the electrons are now integrated out to ob- 
tain an effective action for the order parameter, gen- 
erating a priori infinitely many interaction terms. We 
truncate the effective action, retaining the leading fre- 
quency and momentum dependence of the Gaussian part 
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of the action as well as a constant interaction term, 
since all other terms are irrelevant in the renormaliza- 
tion group sense^'^'^^ (cubic terms are discussed in ap- 
pendix A). For the quadratic part one obtains 5*2 = 

^E.fc*l,fe('5aa'/^ + xL'(fc>i^n))*-a,„,-fc where J is 



the interaction in the spin-spin channel and x^^^, (fc, ia;„) 
the susceptibility in the presence of the finite field B eval- 
uated at J = 0. Calculating these susceptibilities on the 
paramagnetic side of the transition we obtain 



5 = 52[*]+^4[*], (4) 
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Here /3 = l/ksT and w„ = 2nn/l3 is a Matsubara fre- 
quency and fc is measured again from the ordering wave 
vector. The coefficients of and |a;„|cos^ have been 
made unity by an appropriate choice of the bare length 
scale ^0 and temperature/energy scale Tq. In general the 
prefactors of the fc^ and |a;| terms for $2 and ^x/y will be 
different (even after rescaling), we suppress these prefac- 
tors to keep the notations simple as they will not lead 
to any qualitative changes in the results. It is, however, 
essential to keep track of the dynamics of ^x/y, i-S- of 
the ratio of precession and damping terms parametrized 
by an angle 9. For small the dynamics is overdamped, 
while for 9 ~ tt/2 precession dominates. The value of 9 
depends on details of the band-structure and the size of 
the magnetic field with 9 oc B for small magnetic fields. 

As anticipated in (3), the x— and y— direction are 
coupled for 9 > 0. The Gaussian part of the action 
can be diagonalized by introducing the complex field 
<1?^ = ($a; -|- i$y)/-\/2 as above, and we obtain 

S2[$\$1 = J ^l^$^^,,*(2xfc(ia;„)-^)$;^„,, 

+ K^A>''+r. + \u;n\W-^^,-k, (6) 

where 

Xfe(ia;„) = (fc^ -I- r_L -I- |a;„| cos 6* - iuj„ sin^)~^ (7) 

is the propagator of 

As expected from the symmetry arguments given 
above, r± and turn out to be different with > r± 
and rz — r± oc for small B. As r^/^ increases for 
increasing fields {r{B) « r(0) +cB^ for small B), an an- 
tiferromagnetic system sufficiently close to its quantum 
critical point can be tuned to the paramagnetic phase by 
applying a magnetic field (assuming that no first order 
transition is induced). 

When discussing the behavior close to the quantum 
critical point, it is important to note that the mag- 



netic field enters into the calculations both by the in- 
dependence of r as well as through the 5-dependent angle 
9. Close to the quantum critical point tuned by a finite 
magnetic field Be, 9{B) w 9{Bc) can be approximated 
by a constant (as checked below) while it is obviously es- 
sential to keep track of the leading B dependence of the 
control parameter r{B) cc B — Be- 

At this point, it is worthwhile to take a closer look 
at 52 in coordinate and time space for 6 = 7r/2, 
i.e. if Landau damping is absent as it is the case in an 
insulator like TlCuCla* or in a metal with Q > ikp (see 
introduction). The Gaussian part of ^2 is minimized 
for a field that obeys the equation 

idt^^ = H^^, H ={-W^ + ri_). (8) 

This has the form of a Schrodinger equation for a particle 
in a constant potential given by ^ = . If one adds the 
interactions one obtains a non-linear Schrodinger equa- 
tion or Gross-Pitaevskii equation which describes the 
physics of weakly interacting Bosons. In this interpre- 
tation, r takes over the role of the chemical potential. 
The quantum critical point of a field tuned insulating 
antiferromagnet {6 = 7r/2) is therefore in the same uni- 
versality class as the quantum phase transition of a di- 
lute gas of Bosons driven by a chemical potential. The 
non-magnetic phase (r > 0) corresponds to a phase with 
negative chemical potential where no Bosons are present 
in the T ^ limit, while the Bose-Einstein condensed 
phase corresponds to the magnetically ordered phase. 



II. RENORMALIZATION GROUP EQUATIONS 
AND CORRELATION LENGTH 

The physical properties of the eflFective action (4) can 
be analyzed with the help of renormalization group equa- 
tions. As a first step it is useful to perform a simple 
scaling analysis of S{^\. When momenta, frequencies 
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and fields are rescalcd as k' = kb, lu' = Lub^, where z 
is the dynamical critical exponent, and = $6" ''^2'^^ , 
S[^] remains invariant under scaling provided that z = 2. 
The masses r-^'^ and the dimcnsionlcss coupling constant 
u = 5^o/-^o have the scaling dimensions 2 and 4— (d+z), 
respectively. In an aiitiferromagnetic metal, damping as 
well as precession are linear in frequency and the terms 
therefore behave in the same way under scaling. In the 
renormalization group terminology this implies that the 
precession term is an "exactly marginal" perturbation 
with respect to Hertz fix-point {6 = 0,u = 0) which can 
be expected to modify the behavior of the system at the 
quantum critical point. 

The renormalization group equations for the parame- 
ters T, r, and u with corrections to scaling can be ob- 
tained by closely following the procedure employed by 
Millis^'': We introduce a UV-cutoff in the linked cluster 
expansion of the free energy and express changes of that 
cutoff in terms of changes of the parameters of the model. 
The RG equations are as follows: 



dTjb) 
dlogb 
dr^jb) 
dlogb 

dr,{b) 
dlogb 

dujb) 
dlogb 



zT{b), (9) 

2r^{b)+Mb){V2{r±{b),T{b)) 

+ fi{rAb),T{b))), (10) 

2r,{b)+Au{b){f^{r^{b),T{b)) 

+ ifl{r,Xb),r{b))), (11) 

(4-(d + ^))w(6), (12) 



where T is the running temperature and the expressions 
for J2'^' as well as details of the calculation can be found 
in appendix B. Since the scaling dimension for u is neg- 
ative for an antiferromagnetic system in 3 spatial dimen- 
sions, we only consider contributions up to and including 
first order in u. To this order, the scaling law for u re- 
mains unmodified, and 6 remains unrenormalized. The 
parameter obtains, however, finite corrections by higher 
order contributions. 

Equations (9) and (12) are solved trivially. As (6) > 
rx {b) , $2 remains massive at the quantum critical point 
(r^ > for rx =0). In the following wo will concentrate 
on the regime T < , where the influcnc;c of the; parallel 
mode $2 can be absorbed in a redefinition of the bare 
r_L. 

Eq. (10) can be solved for low temperatures in the lim- 
its r_\_/T <C 1 and r_\_/T ^ 1, corresponding to quantum 
critical and (renormalized) Fermi liquid regime, respec- 
tively (see Fig. 1). This provides us with an expression 
for the correlation length We refer to appendix B for 
details of the calculation. In the quantum critical regime 



^x is given by 

^I^(r-± < T) = rx + 16\^7r^/^ C(3/2) uT^'^ cos{9/2), 

(13) 

and in the Fermi liquid regime it has the form 



^l\T « rx) = rx + ^^\TVf 



cost 



(14) 



For all < 7r/2 one obtains the same qualitative behavior 
as in the case of vanishing external magnetic field"'^'*. Only 
in the Fermi liquid regime for 6 = n/2, the correction 
is suppressed as Landau damping is absent in this limit 
and our model is characterized by an energy gap which 
leads to an exponential dependence exp(— r_L/T) of the 
correlation length. 



III. THERMODYNAMIC QUANTITIES 

In this section we calculate the specific heat 7, the 
temperature dependence of the magnetization, the mag- 
netocaloric effect F^, and the susceptibility. The free 
energy can be calculated directly from RG equations fol- 
lowing again Ref. [14]. However, as the quartic coupling 
u is irrelevant, the leading behavior in the paramagnetic 
phase can equivalently be extracted just from the Gaus- 
sian free energy 

T= - F{T = 0)) 



ToV 

'ii (2^ 



cPk f (kij 



coth ( 

2 



- 1 



arctan 



2(r + P)lo cose 
(r + fc2)2 - ^2 



(15) 



measured in units of TqV/^q, and we have set r = r_\_. 

In Eq. (15) and in the results shown below we ig- 
nore contributions from the massive, non-critical mode 
$2 characterized by a finite mass r^. To leading order, 
the corresponding (analytic) corrections to the free en- 
ergy and its derivative are just additive and can be ob- 
tained by replacing r by r^, by setting 6 — Q and by 
dividing the result by a factor 2 (as there are two modes 
perpendicular to B) in all formulas for thermodynamic 
quantities given below. 



A. Specific heat 

We first consider the specific heat coefficient cy/T 
7(T, r) = —d'^T jdT"^ . More precisely, we calculate 7 
7(r,r)-7(T = 0,r = 0): 
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7=4 



(Pk f dx 2x 



(27r)3 



+ r)3 (4(fc2 + r)2((fc2 + + 2rV) cos 61 + 4TV cos(36I)) 4cos61 
((fc2 + r)4 + T4a;4 + 2(A;2 + r)2T2a;2 008(261))^ 



(16) 
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FIG. 2: Scaling function for the specific heat coefficient 
-^7, where a non-critical contribution has been subtracted, 

7 = 7(r, r) — 7(r = 0, r = 0). The function is not completely 
universal but depends on the parameter Q. A crossover from 
7 oc ±r^ for T < r to 7 ~ ±v^ for r < T can be observed, 
where the signs depend on the value of d. 7(T) shows a maxi- 
mum for 7r/6 < 9 < 7r/3 as can be seen more clearly in Fig. 3. 



this differs from the physical specific heat by a 
T— independent (but cutoff-dependent) constant 7c cos ^. 

The integrals can be evaluated exactly in certain limits. 
For 6 — n/2 and low T <^ r the specific heat shows 
thermally activated behavior 



7(6* = 7r/2,T^ 0) 



(2^ T372 



exp 



(17) 



as can be expected from a system with a gapped spec- 
trum. For r <C T and T ^ r, i.e. in the quantum critical 
regime and Fermi liquid regime, respectively, we obtain 
for e < ■it/2 



FIG. 3: Specific heat coefficient as a function of temperature 
for = O.SiT and different values of r. Note that the total 

specific heat coefficient 7 is always positive. The maximum 
for r > is characteristic for systems with n/6 < 6 < it/3. 



negative for 9 < tt/3 and positive for 9 > it/ 3. Also in 
the Fermi liquid regime a sign change can be observed in 
the T"^ /r'l'^ contribution at = 7r/6. 

In Fig. 2 we show the scaling function ^^Jl''^ = fe{T/r) 

obtained from a numerical integration of (16). Due to the 
presence of an exactly marginal perturbation, the scaling 
function is not completely universal but depends on the 
parameter 9. In an intermediate regime, 7r/6 < 9 < tt/3, 
7(T, r) (and the universal scaling function ^^^^=^) shows 

a characteristic maximum as a function of temperature 
as can be read off from the asymptotical results (18) and 
(19). This maximum cannot be seen directly at the quan- 
tum critical point (r = 0) but for any finite r > as long 
as the critical corrections to the specific heat dominate 
the non-critical ones. 



7(r < T) 
7(T < r) 



-V/^cos^ 
6 



60 r3/2 



cos(36'). 



(18) 
(19) 



For ^ = 0, this reproduces well-known results^** (correct- 
ing some factors of 2), and as expected from scaling, ex- 
ponents do not change in the presence of the precession 
term. However, not only the size of the prefactors but 
interestingly also their sign changes when the dynamics 
begins to be dominated by precession rather than damp- 
ing. In the quantum-critical regime the VT correction is 



B. Magnetization, magnetocaloric efFect and 
Griineisen parameter 

As was argued in Ref. [23], the specific heat is not 
the most sensitive thermodynamic quantity close to a 
quantum critical point as it tracks only variations of the 
free energy with respect to temperature (vertical axis 
in Fig. 1) but not with respect to the control param- 
eter B (horizontal axis). It is therefore interesting to 
study also the magnetization M = —dF/dB, the sus- 
ceptibility X = —d^F/dB^, and the T-derivative of M, 
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FIG. 4: Scaling function -^^y- While for the specific heat 

shown in Fig.3 it was necessary to subtract a non-critical con- FIG. 5: |f- as a function of temperature for 9 = O.Stt and 
stant contribution such a background does not exist for different values of r. At the quantum critical point we obtain 

-Vt while ~ -T/VF for T < r. 



OM/dT = -d^F/{dBdT) = -dS/dB. This mixed 
derivative has the advantage that - opposite to specific 
heat coefficient and susceptibility - it vanishes in the 
T ^ limit due to the second law of thermodynamics. 
Therefore it is not necessary to subtract any constant 
non-critical contributions when measuring dM/dT. 

Very interesting are also ratios of the thermodynamic 
derivatives^'^. For a field tuned critical point, one inter- 
esting combination is 



{dM/dT)B 1 {dS/dB)T 



T {dS/dT)B 



1 dT 
T dB 



(20) 

which describes the magnetocaloric effect, i.e. the tem- 
perature change in the sample after an adiabatic change 
of the magnetic field. 

For pressure-tuned quantum phase transitions, where 
d/dB is replaced by d/dp, the quantities related to 
dM/dT, X and Tb are the thermal expansion, the com- 
pressibility (and therefore also to the sound velocity), and 
the Griineisen paramctcr^'^. 

As both control parameter r and 9 depend on B 
one can expect two independent critical contributions to 
dM/dT: 



d'^T dr d^T dO 



dM _ 
'dT ~ ~ dTdr dB 

_ 1 r cfk 



dTdO dB 

TT 



ji2 



dr 



((fc2 + r)4 + uj'^ + 2(fc2 + rfuP- cos(26l) dB 

{IP- ^rfiS}? ^rf -u?)€\^d de 
((fc2 -t- r)2 - a;2)2 + 4(^2 + ^)2^2 ^032 ^5,) ^ 



(21) 



In the limits r, T — > we obtain 



dM 

'Wr^T 



3V27r /3\ /- 



dr 
dB 



dM 

dT T<^r 



IT dr 
127=r'°'^dB 



(22) 
(23) 



As r is a relevant perturbation at the quantum critical 
point while 9 is only marginal, the contributions due to 
the i?-dependence of 9 are subleading and can therefore 
be neglected. 

With r (X B — Be, we therefore find 



dM I- (9 

—— oc — vTcos - 

dM T 

ST T«r ^ ~y/B-B, 



COS 9. 



(24) 
(25) 



Neither {dM / dT)r^T nor {dM / dT)T^r changes sign as 
a function of 9, and indeed, {dM/dT)T,r is a monotonous 
function of T for all values of 9 (see Fig. 4). 

For 9 = n/2 and at low T, the temperature deriva- 
tive of the magnetization also shows thermally activated 
behavior 



dM 

dT 6>=7r/2,T^0 



exp 



(27r)3 VT 

Finally, we evaluate the magnetocaloric effect 

{dM/dT)B 



dr 
dB' 



Tb = - 



T(7 + 7cCOS0)' 



(26) 



(27) 
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which is given by: 



6^/2cos(f)C(i) 



rs(T«r) 



15\/2Tcos(|6l)C(|) - 32\/7V7cCos 61 ' 
1 



2(r - GjcVr) 



(28) 
(29) 



in the limits r,T 0. Due to the non-critical contribu- 
tion 7c, the result for T — > is not fully universaP^. 



C. Susceptibility 

Since r ~ B — Be and T have the same scaling ex- 
ponents, one might expect that the susceptibility x = 
dM/dB = —d'^F/dB'^ and the specific heat coefficient 
7 = —d'^F/dT'^ show very similar behavior. This, how- 
ever, turns out to be not correct. The technical reason 
for this is that the susceptibility in the quantum criti- 
cal regime is a singular function of the (dangerously ir- 
relevant) spin-spin interaction u. Practically, this im- 
plies that a measurement of the susceptibility is comple- 
mentary to other thermodynamic measurements as it is 
highly sensitive to a quantity which can otherwise be de- 
termined only by neutron scattering measurements of the 
correlation length. 

The susceptibility x = —d'^F/dB'^ gets contributions 
both from the B-field dependence of 6 and of r. We only 
consider the leading corrections due to r (see discussion 
above) and evaluate the quantity x{r,T) = x{'r,T) — 
X{r = 0, T = 0) with 



+ Q{-Lu)lm 



Hb (w)Im 



1 



{k^ + r + iw cos — wsm6)^ 



(A;2 + ito cos 6 — LOsinO)^. 



(30) 



Most interesting is the quantum critical regime, where 
the leading correction to the susceptibility takes the form 



X{r « T) 




47r2 



VTcos 




(31) 



Note that this expression formally diverges for r 
0. This implies that we have to take into account the 
interaction effects discussed in Section II and we have to 
replace the control parameter r by ^J^{T) r + uT^/"^ 
given by Eq. (13). One therefore finds 



X oc for r < T < {r/uf^, (32) 



but 



X oc for T > (r/uf/^. 



(33) 
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FIG. 6: Susceptibility as a function of temperature for u = 
0.2, Q = 0.01, and different values of r. Curves for other values 
of 9 look essentially identical and are not shown. Note that the 
susceptibility is much more sensitive to small deviations from 
the quantum critical point than the specific heat coefficient 
or dM/dT (c.f. Fig. 3, 5, where much larger values for r have 
been used). The T^^* cusp at the quantum critical point 
[Eq. 33] is rapidly washed out by tiny deviations from the 
critical magnetic field and replaced by the linear dependence 
of Eq. (32). 



Scaling is violated, i.e. the susceptibility is no longer of 
the form x(r, T) = f${T/r)/ V^, as the dangerously ir- 
relevant coupling u determines x in this regime. It is 
interesting to trace back the origin of the l/\/r contri- 
bution in (31). It arises from the w„ = mode of the 
Gaussian theory (6) which leads to a contribution of the 

form TJ2k (^w+r) *° susceptibility. Note that the 
static LU = contribution does not depend on the dynam- 
ics, i.e. it does not depend on 9. However, the correlation 
length ^_l(T) given in Eq. (13) does depend smoothly on 
6 which leads to a slight ^-dependence of 

In the Fermi liquid regime, the susceptibility is given 

by 



dr 
dB 



x(r«r)=(-^yFcose+l^cos^) ( 

(34) 

while for 6 = ^/2 and low T ^ r it shows thermally 
activated behavior 



X(0 = 7r/2,T^G) 



(27r) 



exp 



dr 
OB 



(35) 

The rapid crossovers between the T^, T and T^'* 
regimes are shown in Fig. 6. 



IV. SCATTERING RATE 

Due to energy and momentum conservation, the scat- 
tering of electrons from spin-fluctuations is most efS- 
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cient close to "hot lines" on the Fermi surface, where 
Ek = Ek±Q = 0, where Q is the ordering vector of the 
antiferromagnet. In order to determine how spin preces- 
sion modifies the results we calculate the scattering rate 
as a function of 9. We will neglect all orbital effects of 
the magnetic field and will not try to calculate the con- 
ductivity and the Hall effect. For an extensive discussion 
of orbital effects and magnetotransport in nearly antifer- 
romagnetic metals see Ref. [22] which does, however, not 
consider effects of spin precession. 

In second order perturbation theory the lifetime of the 
spin-up electron scattering from fluctuations of at 
T = is given by^^ 

4 = E / ''(}uImxk-k'{uj)6[uj-{E+-E-,)], (36) 

I 



where is a coupling constant, i?^ and Vp the en- 
ergy and velocity of spin up/down electrons and x is the 
spin fluctuation spectrum of Eq. (7), 



Xq{i^) = : —■ — ^—5 

LOq + r + ICO COS — cosmO 



with LOq = {q±Q)^ /qq. We split the momentum integra- 
tion in an integral over the Fermi surface and an energy 
integration J (fk' = JJ dk' /vp J dE^, and integrate first 
over E^, , then over w to obtain 



f(2^)' 

2 2 

9s% 



dk' [ cos 6 In 



:E1 



cos ( 



(a;fe_fe'+r)2+2£+sin0 + (i;+)2 



(wfe.fe' -h rf 



+ sin Q arctan 



-E^cose 
.k' +r + E+ sine _ 



-0 



(38) 
(39) 



where Sk = r + (Sk/qQ)"^ and Sk is the distance of fc + Q 
from the Fermi surface or, approximately, the distance of 
k from hot lines on the Fermi surface. Analogously we 
obtain for spin-down electrons 

4 - -rjf^E^^^l^^ose, I -9], (40) 

t;J(27r)2 [2^^ 2 J 

where the indices + and — have been exchanged w.r.t. 
(38). The scattering rate is strongly dependent on the 
distance from the hot lines; at the quantum critical 
point and for Sk/qo fa the scattering rate is linear in 
the quasiparticle energy. Far away from the hot lines 
and the quantum critical point the usual scattering rate 

l/r^"'' oc E^'~ is recovered^^. The main result of this 
section is that the spin-precession term does not lead to 
a qualitative change in the scattering rate. 



V. DISCUSSION 

In this paper, we have discussed the field-induced 
quantum phase transition of a clean, three-dimensional 
antiferromagnetic metal, restricting our attention to the 
non-magnetic side of the phase diagram. The main ques- 
tion was how the interplay of precession of the spins in 
the presence of a finite magnetic field and Landau damp- 
ing modifies the quantum critical behavior. 

One main qualitative result of our analysis is that the 
critical behavior is not completely universal as it depends 
on a continuous variable 9 which parametrizes the ratio 



of precession and damping terms in the effective action. 
While critical exponents do not depend on 9, this param- 
eter strongly changes the scaling functions and even the 
sign of leading corrections e.g. to the specific heat. 

A requirement for the validity of our analysis is that 
the two modes ^x,y perpendicular to the magnetic field 
are characterized by the same mass. In the presence of 
sizable spin-orbit couplings, this will only be the case, if 
the crystal has a sufficiently high symmetry and if fur- 
thermore the external magnetic field is applied along a 
symmetry direction of a crystal. 

Presently, we are not aware of any experiments which 
show for example the maximum in the T dependence of 
the specific heat coefficient which we predict for n/Q < 
9 < 7r/3. Note that in systems like CeCu5.2Ago.8^ or 
CeCu5.8Auo.2^ strong anisotropics prohibit the preces- 
sion of the spin, i.e. 9 — i). Under what conditions can 
large values of 9 be expected? Obviously, large uniform 
magnetizations are required. In heavy Fermion systems 
with Kondo temperatures of the order of a few Kelvin, 
one can introduce strong magnetic polarizations with 
moderate external fields and it should therefore be pos- 
sible to induce sizable values of 9. A different class of 
systems which might be of interest in this context are 
ferrimagnetic materials. If it is possible to suppress only 
the staggered component of the magnetization in such 
systems either by external fields, pressure or doping, the 
critical theory within the Hertz approach will be charac- 
terized by a finite (and again sizable) 9. 

As long as the ordering vector Q of the antiferromagnet 
can connect the spin-up and spin-down Fermi surfaces 
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control parameter r 



FIG. 7: Qualitative behavior of specific lieat coefficiont 7 and 
magnetic susceptibility x in various regimes of the phase dia- 
gram. 

{Q < kp + kp or, more precisely, E'^ = EJ^^q = for 
a line of momenta fe), Landau damping is present and 
6 < 7r/2. In contrast, one finds ^ = 7r/2 in all systems 
where no such connection exists {Q > kp+kp). However, 
the transition from < 7r/2 to = 7r/2 is not expected 
to be smooth, as the interactions of the spin-fluctuations 
diverge and become relevant^^'^^ at the point where Q = 
kp -]r kp . 

According to our analysis, the susceptibility dM/dB is 
a particularly interesting experimental quantity to study 
close to a field-driven quantum critical point. First of all, 
it is expected to be much more sensitive to small devia- 
tions from criticality compared to other thermodynamic 
quantities (see Fig. 6). Second, it allows to measure the 
correlation length, a quantity which cannot be extracted 
from other thermodynamic quantities, as for B = we 
obtain from (31) 

Third, it strongly violates the T/ (B — Be) scaling. This 
deviation from scaling for x can be used to show that the 
relevant critical theory is above its upper critical dimen- 
sion, a central question for the interpretation of quantum 
criticality in systems like CcCue-jjAua: or YbRh2Si2^'''. 
All these three statements arc actually independent of 
the value of 9: they apply equally for a dynamics which 
is overdamped, ^ -C 1, or for a BEC system like TlCuCla 
with 9 = 7r/2. Note that in pressure-tuned quantum crit- 
ical points the compressibility k (and therefore the sound 
velocity^^') plays the same role as % for field tuned quan- 
tum phase transitions. An overview of the qualitative 
T dependence of the specific heat and susceptibility is 
shown in Fig. 7. 

Field tuned quantum-phase transitions in metals allow 
to study quantum critical behavior with a tuning param- 



eter which can easily be controlled and with a conjugate 
field - the uniform magnetization which can directly be 
measured. They are therefore especially well suited to an- 
swer some of the central questions in the field of quantum 
critical metals, for example, whether or not such systems 
can be described in terms of simple spin-fluctuation the- 
ories as have been used in this paper. 
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APPENDIX A: CUBIC TERMS IN THE 
EFFECTIVE ACTION 

In this appendix we briefly discuss whether cubic terms 

arc present in the low-energy effective Lagrangian. As 
$ carries the momentum Q, the presence of such terms 
is forbidden by momentum conservation in most systems 
with the exception of magnetic structures (e.g. BCC lat- 
tices) where the sum of three ordering vectors adds up to 
0. If such a system has Ising symmetry, then a cubic term 
does exist and the magnetic field driven transition will be 
first order. However, for xy symmetry perpendicular to 
the magnetic field (the case mostly discussed in this pa- 
per), a rotationally invariant cubic term of the form 
does not exist. While terms like B^z\^±\'^ are allowed 
by symmetry, they lead effectively only to a rcnormaliza- 
tion of the interaction as <l>z remains massive. We 

therefore neglect such terms. 



APPENDIX B: DERIVATION OF 
RG-EQUATIONS 

Following Millis' treatment^^'^^, we perform the renor- 
malization group analysis on the free energy after hav- 
ing converted all Matsubara sums to integrals. Although 
we restrict our calculations to systems in three spatial 
dimensions and with a dynamical critical exponent of 
z = 2, we nonetheless keep the variables d and z in the 
calculation in order to make the origin of certain factors 
more transparent. 

The free energy can be obtained via a linked cluster 
expansion in the coupling constant u. The scaling di- 
mension of u is 4: — {d + z), which is negative for an an- 
tiferromagnetic system in 3 spatial dimensions. To first 
order in ?i, only the diagram in Fig. 8 contributes to the 
free energy. Up to first order in u, the free energy is 
therefore given by 

J- = J-G + u[{I^ + J^. +I^.f+2{ll + J2. + )] , (Bl) 
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+ 




FIG. 8: Diagram contributing to the free energy in 0(u). 
The different contractions of the internal indices of the fields 
involved have been made explicit on the right hand side, where 
the dashed line represents the quartic interaction u. 



where 



A r 

Tg = / 7 TT / — coth — X 

2 J (27r)3 J n \2 ' 







arctan 



2(r + P)ujcos( 
(r + - ^2 



(B2) 



is the Gaussian free energy measured in units of ToF/^q 
with the cutoffs A and F, is given by 



A r 

d^k f dio , f (iuo , 
,„ ,„ / — coth -— X 
(2^)3 J n V 2 





(r + fc2)2 _ 2(r + B)uj sin 6* + 



(B3) 



and I^* , I^z are defined analogously. 

As a next step, we separate out of the momentum and 
frequency integrals in the expressions on the right hand 
side of (Bl) the regions given by {A > fc > A/6, F > > 
0} and {A > > 0, F > w > V/h'^). Using that 



7, + /,=2|^2.I,=2^. (B4) 



the change in J- upon such a variation of the cutoff can 
be expressed as a change of r± and Tz, and this leads to 
the equations 



9 log 6 
9 log 6 



= 2r^[h)+Au{h){2f^{r^{b),T{h)) 

+ /|(r,(fe),r(fe))), 
= 2rz(6)+4u(6)(/2^(r^(&),T(6)) 



(B5) 



+ 3/|(rz(6),r(6))) (B6) 
for the running masses r±{h), rz{b), where and /| are 



given by 



/-K,T)=ir3A3/^coth(^ 



2cj((A2 + r^)2 +Lu'^) cost 



((A2 + ri)2 + cj2)2 _ 4(A2 + r_L)2a;2 sin' 

A 



zT 



d^k , //3F\ 
- — -- coth — X 
(27r)3 V 2 ; 

2F((fc2 + r_L)2 +r2)cos( 



dio 

TT 



r2)2 - 

coth 



4(fc 
/3w 



r_L)2F2 sin' 6* 



UJ 



zT 
H 

TT 



d^fc , //3F 
- — coth — 
(27r)3 V 2 



(r^ + A2)2 +tj2 

F 

(rz+fc2)2+r2- 
(B7) 



In the following we assume that the system is close to 
the quantum critical point at temperatures much smaller 
than r^. In this case, fii'^'z,'^) can be set to zero, and 
the renormalization group flow of r± is determined by 
{r± , T) only. There are two contributions to , one 
from the renormalization due to the separated momen- 
tum shell, where momentum is set on shell fc = A, and 
one from the renormalization due to the frequency shell 
with UJ ^ T. For subsequent calculations we note that 



f2ir±,r)-f^{r^,0) 



djjj 

TT 



coth 



(3uj 



2w((A2 +rj_f +w2)cosf 



((A2 + r^)2 + tj2)2 _ 4(^2 r±)^uj^ sin' 



+ 0(e-i^/^), 
(B8) 



in other words the contribution of the frequency shell 
renormalizes zero temperature properties only and is ex- 
ponentially suppressed at finite temperatures. 

In order to obtain an expression for the correlation 
length, we first substitute r±{b) = R±{b)b^ to eliminate 
the naive scaling and then formally integrate equation 
(B5): 



Inb 

8 J dxe~ 




■2a; „ 



We then perform an expansion in temperature 

R^{b) ^ ^Ab) + Rrih) + 5RL{b), 
where three terms contribute. 



(B9) 
(BIO) 
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The first term A±{b) is the running mass at zero tem- and SR±{b) originates from the temperature dependence 



perature, 



In 6 



Ai(6) = ri+S j dxe-2-«(e^)/2^(Ai(e^)e2-,0); 

(Bll) 

the integrand can now be expanded in Aj_ which leads 
to the following expression 



In 6 



Ai(6) ~ ro^ + 8/2-^(0,0) j dxe-^^'uie') 



^0 ~ 



(B12) 



i.e. this defines the parameter which characterizes the 
distance from the critical point. 

The other two terms in (BIO) are of first order in tem- 
perature: One contribution is due to an explicit depen- 
dence of f2 on the running temperature, 

In 5 

R^{b) = 8 ydxe-2^u(e^)(/2-^(i?_L(e^)e2^,Te^^) 



of the running mass, 



In 6 



5R^{b) = 8 y"da;e-2^u(e^)(/2-^(i?_L(e^)e2^,0) 



(B14) 



This term is of the order of u and will be neglected from 
now on. 



(B13) by 
I 



The inverse square of the correlation length is given 



Inb 



= lim {A(6) + R^{b)} = r±+ \im8 [ dx e-'^-'=u{e'')K^h^ x 



IT 



coth 



V2Te^^/ 



2u;((A2 + i?i(e^)e2^)2 + tu^) cos 



(B15) 
(B16) 



((A2 + i?x(e^)e2^)2 + a;2)2 _ 4(^2 + i?_L(e^)e2^)2a;2 sin^ 6 

U 

/■OO POO 1 

+ leA-^+^-^ifdr^/^ / (ia;u(e"AT-i/^)e("-")^ / — (cothz; - 1) 

Jhi{^) Jo 

4A^e^^t;((A2 + ii_L(e^ATV^)e2^A2T-2/^)2 + AA'^^e^^^v'^) cosO 

I 



where the transformations e^ = e^A ^T^/^ and v = 
uj/2Te^^ have been introduced, and u(e^AT~^/^) = 
uo(e^AT-i/^)4-('^+^). Expression (B16) for the correla- 
tion length can now be evaluated in the quantum critical 
and Fermi liquid regime. 

In the quantum critical regime we can neglect the de- 
pendence of the integrand of (B16) on R± and extend the 
lower limit of the a;-integral to —00. Using the following 
integral 



f 

Jo 



rfC ^4t7 = 2nr(n)C(n), n = 0, 1, 2, 
smh 4 



(B17) 



we obtain 



r_L + 16 



d + z-2 

xuT ^ cos 



^cos(i^^) 

rf- 2 



(B18) 



In the Fermi liquid regime and for low temperatures, 

we can replace the running mass i?_L in (B16) by the 
control parameter r^. It is convenient at this point to 
introduce yet another variable transformation of the form 
e^^ = r xT^"^! . To lowest order we can then neglect 

— z/2 

the term Trj_ in the integrand. Furthermore, we can 
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extend the lower limit of the a;- integral to — oo, thereby in the Fermi liquid regime, 
inducing an error of order 0{rf^ / K)'^~'^'^^ , and obtain 



Cl'=r^ + 16^ . . KauT\'~- cosO (B19) 

1^ sm (^^vrj 
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